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Abstract 

The energy dependence of the neutrino-iron and antineutrino-iron inclusive charged-current cross 
sections and their ratio have been measured using a high-statistics sample with the MINOS Near 
Detector exposed to the NuMI beam from the Main Injector at Fermilab. Neutrino and antineutrino 
fluxes were determined using a low hadronic energy subsample of charged-current events. We report 
measurements of u-Fe (u-Fe) cross section in the energy range 3-50 GeV (5-50 GeV) with precision 
of 2-8% (3-9%) and their ratio which is measured with precision 2-8%. The data set spans the 
region from low energy, where accurate measurements are sparse, up to the high-energy scaling 
region where the cross section is well understood. 
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I. INTRODUCTION 



Neutrino-nucleon and antineutrino-nucleon charged- current [y^N CC and z/^N CC) in- 
clusive cross sections above 30 GeV have been determined by several experiments jj-y] with 
a combined precision of 2% [4]. The measured cross sections at these energies have a linear 
dependence on energy, which agrees well with the prediction of the Quark Parton Model 

(QPM) Q. 

At lower energies, the cross section is both less well measured and difficult to model 
due to overlapping contributions from quasi-elastic processes (u^ + n — > /i~ + p), resonance 
excitation followed by subsequent decay, and the onset of deeply-inelastic scattering (DIS). 
This energy range is of particular interest to ongoing and future neutrino oscillation searches 
in MINOS, NOz/A jfj], and T2K Most cross section measurements in the E v <30GeV 



range 



14j have uncertainties of the order of 10%. Recently, NOMAD [15j has measured 
the cross section down to 2.5 GeV with a precision of better than 4%. However, this result 
relies on a particle production model tuned to data [16| to predict the neutrino flux. In this 
paper we present a measurement of the v^N CC cross section with a precision from 2-8%, 
covering the 3-50 GeV energy range using the MINOS Near Detector. Our analysis uses a 
low hadronic energy subsample to determine the flux shape [l?j, 18]. 

Antineutrino-nucleon charged-current cross sections in the E v <30 GeV range suffer from 
the same complications listed above and tend to be even less well measured. Several ex- 



periments reported results 



11 



13 



19 



20]; however data coverage in energy was sparse and 



these measurements typically have larger than 10% uncertainty. Our measurement has 
higher precision, with uncertainties which range from 3-9%. 

The z/ M N CC to z/ M N CC cross section ratio, r = a u /a u , has been measured with a 



measurement 



combined precision of better than 1% at high energies [17, |2l| but only one dedicated 



22( | has been performed in the E u <30 GeV range. Gargamelle [22J reports 



measurements of r from 1-10 GeV with precision of about 20%. Our result substantially 
adds both coverage and precision to the determination of r. The ratio is more precisely 
determined than either cross section measured separately due to a partial cancellation of 
most systematic effects and a cancellation of the normalization uncertainty. 

The results in this paper can be used to tune and improve neutrino interaction generator 
models 23j, |24j. For example, neutrino scattering data are required for the modeling of 
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251 ] . Also, the cross section ratio r is 



the axial vector contribution to the cross section 
particularly sensitive to the modeling of xF%, the parity violating structure function, which 
enters into the numerator and denominator with opposite sign, and to the antiquark content 
of the nucleon, which contributes differently to neutrino and antineutrino scattering. In 
addition, at 5 GeV, about 70% of our event sample has negative 4-momentum transfer 
squared, Q 2 , of less than 1.5 GeV 2 . This large, low-Q 2 sample provides model sensitivity to 
the low-Q 2 QCD contributions (higher order QCD, higher-twist, and target mass corrections) 
which are difficult to calculate. 



Overview of the Analysis 

The Ufj, CC and CC total cross sections as a function of incoming neutrino energy E 
are determined from the inclusive charged- current interaction rate and the incident neutrino 
flux. A sample of CC events ("cross section sample") is selected and a subsample of these 
events with low hadronic energy ( "flux sample" ) is defined. A Monte Carlo simulation which 
includes detailed detector geometry and response is used to correct the flux and cross section 
samples for detector acceptance and smearing effects. 

Neutrino and antineutrino differential cross sections, ^j— , approach the same constant 
value, independent of energy, in the limit of low-z/, where v is the energy transfered to the 
hadronic system. A method which exploits this feature is used to determine the energy 
dependence of the flux from the flux sample, which is then normalized using the world 
average cross section value measured above 30 GeV. To accomplish this we make use of the 
full range of our data sample, which overlaps with the high energy measurements in the 
30-50 GeV region. This "low-^" method has been used previously at high energies jl?] . rsl ] 
and here it is adapted to the E <30 GeV range. 

The neutrino beam, detector and the Monte Carlo simulation of the experiment are 
described in Sec. [Til Sec. II III describes the event sample selection and the methods for 
extracting the flux and the cross section. A discussion of systematic uncertainties and 
results are given in Sec. [TV] and [V] respectively. 
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II. BEAMLINE AND DETECTOR 



MINOS is a two-detector, long baseline neutrino oscillation experiment using the NuMI 
(Neutrinos at Main Injector) neutrino beam at Fermilab. The oscillation parameters are 
measured [26, 0] by comparing the energy spectra at the Near Detector located at Fer- 
milab and the Far Detector located 734 km away in the Soudan Mine in northern Minnesota. 
In this section we describe the neutrino beam, the Near Detector, and the Monte Carlo sim- 
ulation. More detailed descriptions of the beamline and the MINOS detectors are given 



elsewhere 



28]. 



A. Neutrino Beam 



The NuMI neutrino beam is produced from 120 GeV protons extracted in a 10 /is spill 
from the Main Injector which impinge on a graphite target, with a typical intensity for 
the data presented here of 2.2 x 10 13 protons on target (PoT) per spill. Charged particles 
produced in the target, mainly pions and kaons, are focused by a pair of toriodal magnets 
called horns into a 675 m long decay volume where the mesons decay to muons and neutrinos. 
The decay region is followed by a hadron absorber where remaining mesons and protons are 
stopped. The neutrino beam then traverses 240 m of unexcavated rock before reaching the 
Near Detector located 1.04 km from the target. 

Data for this analysis were collected in "low energy" beam mode in which the downstream 
end of the target is placed 10 cm from the neck of the first focusing horn and the current in 
the horns is 185 kA, with the polarity set to focus positively charged mesons. The Monte 
Carlo simulation predicts the composition of the event sample to be 92.9% 2^, 5.8% u^, and 
1.3% v e + i> e . Fig. CQshows the simulated flux spectrum of the and in the beam. The 
component of the beam, which results primarily from focused ir + and K + , peaks between 3 
and 4 GeV with a long tail. The component arises mainly from low transverse momentum 
7r~ and K~ traveling through the neck of both horns, where they undergo little defocusing. 
This results in a spectrum with no focusing peak and greater mean energy. 
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FIG. 1: The muon neutrino and antineutrino flux at the center of the Near Detector as calculated 
by the NuMI beam simulation. 

B. Near Detector 

The Near Detector is a tracking calorimeter composed of planes of magnetized iron and 
plastic scintillator. A toroidal magnetic field with an average strength of 1.3 T provides 
a measure of muon momentum from curvature and is used to distinguish and CC 
interactions based on the charge sign of the final state muon. In normal operational mode 
the field is set to focus negative muons. 

The Near Detector, illustrated in Fig. [2j consists of 282 steel plates, 2.54 cm thick, of 
which 152 are instrumented with 1 cm thick scintillator planes. The scintillator planes are 
made of 4.1cm wide strips oriented ±45° with respect to the vertical and alternating ±90° 
in successive planes. The strips are read out with wavelength shifting fibers connected to 
multi-anode photo- multiplier tubes (PMT). Every fifth plane throughout the detector is 
fully instrumented with a scintillator layer. In the upstream calorimeter region, comprising 
the first 120 planes, each of the four intervening planes has partial scintillator coverage. The 
calorimeter region is used to measure energy deposited by neutrino-induced hadronic show- 
ers. Event verticies are required to be within a fiducial volume contained in the calorimeter. 
The downstream 162 planes of the detector form the muon spectrometer. 

In the low energy NuMI beam configuration, the typical interaction rate in the Near 
Detector is about 16 events in a 10 /is spill. Events are separated using timing and spatial 
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FIG. 2: Left: top view of the Near Detector, showing the calorimeter and muon spectrometer. The 
drawing is not to scale. Right: transverse view of a Near Detector plane. The shaded area shows 
a partially instrumented active scintillator plane and the dashed line within shows the boundary 
of the fiducial region. The dotted line shows the outline of a fully instrumented scintillator plane. 

information. The events accepted for this analysis were from interactions occurring during a 
13 /is long gate synchronized to the beam spill. The readout electronics continuously digitize 
the PMT signals in 19ns samples without deadtime throughout the spill. In between beam 
spills, cosmic ray muon data are recorded with less than 1% deadtime. 

The detector is calibrated in several steps that convert the raw PMT signal to deposited 
energy 



281 ]. The non-linearity of the electronics is measured with charge injection; relative 
PMT gains are measured with an in-situ light injection system; variations in the light output 
between scintillator strips and along the strips are corrected with cosmic ray muons and a 
radioactive source scanner. Cosmic ray muons which stop in the detector are used to calibrate 
the measured signal to energy lost by muons passing through the scintillator strips. The 
detector simulation is tuned to emulate the actual detector response at all stages in the 
calibration chain. 



C. Beam and Detector Monte Carlo Simulation 

A Monte Carlo simulation is used to model the production of the neutrino beam, interac- 
tion of neutrinos in and around the detector, and the detector response, which is simulated 



using GEANT3 



291 ] . The beam model includes a simulation of secondary hadron production 
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from proton interactions 30] and the propagation of these hadrons. Their reinteraction 
and decay products are also tracked through the target, magnetic horns, and decay region. 
This simulation produces an initial estimate of the flux, which is later replaced by the flux 
extracted using the method described below. 

Neutrino interactions in the detector are simulated using the NEUGEN3 [23J event gen- 
erator. The simulation of quasi-elastic interactions, which dominate at low energies, is 



based on the Llewellyn-Smith 



31| model, while intermediate-energy resonance interactions 



are simulated according to the Rein-Sehgal model [32j, |33j. Both models assume a dipole 
parametrization of the axial part of the cross section that depends on the axial mass param- 
eters M A (QEL) and M A (RES), taken to be 0.99±0.15 and 1.12±0.17GeV, respectively. A 
transition is made between resonance production and the DIS model by phasing out the 
former and phasing in the latter over the hadronic invariant mass range, 1.7< W <2.0 GeV. 
The sum of the resonance and DIS contributions are constrained to match total cross section 
data. 

DIS interactions, which dominate at high energy, are based on an effective leading order 
model by Bodek and Yang [34j. The Bjorken scaling variable x is replaced by an effective 
scaling variable that depends on two parameters Ay lt and B^t, where A^ t accounts for target 
mass effects and higher-twist terms. Bht depends on the transverse momentum of the initial 



34| and gives the parameters 



state quark. The model is fit to charged lepton scattering data 
Aht and Bht and correction factors {C v \ u , C V 2 U , C v \d, C V 2d, Csid and C s i u ) for valence and sea 
up and down quark parton distribution functions. The uncertainties on these parameters 
were not readily available so a study was performed to estimate them and their effect on 
this cross section measurement (see Sec. ITV) . 

The cross section in the transition region from resonance to DIS is expressed as a sum 
of a pure- resonance cross section and a non- resonance contribution from DIS. The sum is 
tuned to describe low multiplicity final state data in this region 23I. For DIS interactions, 
the final state hadronic system is modeled with KNO scaling |35|], which transitions to 
pythia/jetset 36| at hadronic invariant mass W=3 GeV. The total neutrino cross section 
is tuned by a scale factor so that the cross section at 100 GeV matches the world average of 
measurements. 

The dynamics of hadron formation in the target nucleus and reinteraction of hadrons 
after formation modify the visible hadronic shower energy. These effects are simulated using 
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a cascade Monte Carlo anchored to 7rN, pN and 7rFe and pFe scattering; data and validated 
against neutrino-deuterium and neutrino- neon scattering data 37], [38|. A treatment of 
hadron formation time is included [39]. 



III. ANALYSIS 

The CC total cross sections are measured from the inclusive CC scattering rate, Fq-q (E), 
and the incident neutrino flux, <& V ( V \E). A sample of CC events, Nqq (E), is selected 
and then corrected for acceptance and backgrounds to determine T V qq(E). A flux sample, 
F U ^(E), consisting of the subset of Nqq (E) with low u, (in the lab frame v = Ehad, 
the energy measured at the hadronic vertex), is also defined and corrected for acceptance, 
backgrounds, and for a small energy dependence using our Monte Carlo model to yield 
& V ( U )(E). The event reconstruction and selection of these samples to form the cross section 
are described in this section. 

The data used in this analysis were collected between June 2005 and April 2007 and 
correspond to an exposure of 2.45 xlO 20 PoT. The MC sample is almost double the data, 
corresponding to 4.4xl0 20 PoT. 



A. Event Reconstruction 



Neutrino events are identified using the timing and spatial pattern of energy deposited in 
the scintillator strips. Muon tracks are recognized as a string of hit strips typically spanning 
more than 10 steel plates. For muons that stop in the detector the energy is computed from 



range according to the energy loss tables of Groom, et al. [40j . A systematic uncertainty of 
2% is assigned to the energy measured from range, arising from uncertainties in the range 
tables, the variation in material composition and the accuracy of our track length recon- 
struction. The momentum of muons exiting the detector are measured using the curvature 
of their trajectory in the detector's magnetic field. A 4% systematic uncertainty is assigned 
to our knowledge of the absolute muon momentum measurement from curvature. This is 
assessed by comparing the energy measured with curvature to the independent measurement 
from range using tracks that stop in the detector, and by folding in underlying uncertain- 



ties in the detector's magnetic field 



4l| . The resolution for muon momentum measured 
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from range is 5% while that measured from curvature has non-Gaussian tails and width of 
approximately 10%. 

The vertex of a neutrino interaction is taken to be at the start of a reconstructed track. 
Hit strips near the vertex which are not included in the track are identified as coming from 
hadrons produced in the interaction. Their summed signal is converted to energy using a 
lookup table derived from simulated showers to form the hadronic shower energy, Ehad- The 
response of our detector to single hadrons was measured in an exposure of a smaller version 
of the detector to a test beam [42j. The measured test beam detector response was used 



to tune our simulations. The absolute ener gy s cale of the detector's response to hadronic 
particles is modeled to an accuracy of 5.6% 27|, |43), which we take as the hadronic energy 



scale uncertainty in the cross section measurement (see Sec. II Vj) . 



B. CC Event Selection 

The inclusive charged-current sample Nqq (E) is selected using the following criteria: 

1. Fiducial volume - Selected events have a vertex position along the detector axis be- 
tween 0.5 and 4.0 m, measured from the upstream face of the detector. In the plane 
transverse to the detector axis, the vertex is required to be more than 0.5 m from the 
edge of an active scintillator plane and outside of a 0.8 m radius centered at the coil 
hole. The outline of the fiducial region is shown in Fig. |2J 

2. Coil hole - The coil hole is uninstrumented and variations in the material composition 
and magnetic field are somewhat larger in the region around it. To reduce the effect 
of these uncertainties, events with tracks that spend a significant fraction of their 
path-length near the hole are removed from the event sample. A minimum of 95% of 
hit strips in the event are required to be further than 0.3 m from the center at closest 
approach (see Fig. [2]). 

3. Track Energy - The energy of the muon must be greater than 1.5 GeV. This require- 
ment rejects neutral-current (NC) background events, which populate the low energy 
region, and short, poorly reconstructed tracks. 

4. Track Quality - The track fitting procedure yields a measurement of the muon mo- 
mentum with an associated uncertainty. The track fit is required to be convergent and 
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have an uncertainty of less than 30%. In addition, we require the track's longitudinal 
start positions in each view to be less than six planes apart. 

5. Neutrino Energy - The reconstructed neutrino energy, E, which is the sum of the track 
and shower energies, is required to be greater than 3 GeV (5 GeV) for the neutrino 
(antineutrino) sample and less than 50 GeV. The minimum energy requirements are 
imposed to minimize the overlap of the inclusive CC sample and the flux sample, 
which is substantial below these values. Above the maximum energy cut, resolution 
of the track momentum measurement from curvature degrades as the tracks become 
straighter. 

The event sample is divided into two categories depending on whether the track stops in 
or exits the detector. For exiting events, the muon leaves the detector through the back or 
side, or passes into the uninstrumented coil hole region. The stopping and exiting samples 
are further differentiated based on whether they end in the upstream or downstream region 
(see Fig. |2]) because of the difference in sampling in the two regions. 

The CC sample is selected by requiring the sign of the track curvature measurement 
to be positive. This sample has a higher fractional contamination from wrong-sign events 
(misidentified yT tracks) due to the much larger component of the beam. The following 
additional requirements are imposed to the CC sample to reduce this contamination: 

1. Bend away from coil - We require that the track bend away from the magnet coil 
hole to reject positive charge track candidates whose curvature is mismeasured by the 
tracker. An angle is defined in the transverse plane by forming a straight line from 
the extrapolated track end point in absence of a magnetic field to the observed track 
end point, and a line from the magnet coil hole center to the observed interaction 
point {4]]]. For a particle bending toward the coil, this angle will be near it radians, 
while for a defocused y + it will be near or 2ir radians. For the antineutrino sample, 
we select a value for this angle less than 1.04 radians or greater than 5.24 radians. 

2. Number of hit planes - We keep events in which the difference in the number of hit 
planes along the track between the two views is less than five. Events which are 
rejected by this cut usually enter the uninstrumented region in one view yielding an 
unreliable determination of the charge sign. 
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Selection Criterion 


Track Charge < 


Track Charge > 




(%removed) 


(%removed) 


Track vertex in fiducial volume 


3608572 


841986 


>1.5GeV 


2571917 (28.7%) 


344110 (59.1%) 


Track Quality Cut 


2351328 ( 8.6%) 


282657 (17.8%) 


3 < E v <50GeV 


1941019 (17.5%) 




(5 < E v <50GeV for Track Charge>0) 




235024 (16.9%) 


Additional cuts 




159880 (32%) 



TABLE I: Effect of the selection criteria on the negative (left) and positive (right) charge recon- 
structed track samples. Each row shows the number of events remaining after each successive cut. 
The numbers in parentheses show the percentage of events removed by each cut compared with 
the previous row. 

3. Downstream exiting tracks - Only events with tracks that exit the detector in the 
downstream region are used for the antineutrino analysis. The rejected samples have 
high contamination from NC and misidentified CC (wrong-sign) events. 

Table [I] shows the effect of the selection criteria on the neutrino and the antineutrino 
reconstructed samples. The minimum track energy cut has the largest effect, resulting in 
an approximately 30% loss in the sample and 60% in the sample. The cut removes 
primarily NC events, which arise from both neutrinos and antineutrinos, and therefore affects 
the smaller sample more. The track quality cut also has a larger effect on the sample. A 
large fraction of the CC tracks whose charge has been mismeasured due to poor curvature 
determination are removed from the sample by this cut. After all selections have been 
applied, the inclusive event sample Nqq\e) consists of 1.94 x 10 6 and 1.59 x 10 5 
events. 

The CC sample is organized into energy bins and corrected for detector acceptance 
Aqq (E) and backgrounds Bqq (E) to obtain the CC scattering rate, T$$>(E) = (N^P(E)- 
Bqq (E)) I 'Aqq (E) , where A'qq(E), shown in Fig. [31 represents the number of Monte Carlo 
events reconstructed in a given bin divided by the number generated in that bin. The de- 
crease in acceptance at low energy is due to the minimum muon energy requirement. For 
neutrinos, the shape below 10 GeV is determined by the geometry of the detector and over- 
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FIG. 4: Neutral-current and wrong-sign backgrounds in the neutrino (left) and antineutrino (right) 
selected charged-current sample as calculated from the Monte Carlo event sample. The error bars 
show the statistical uncertainty only. 

lap of the stopping and exiting samples, which have different resolutions. The contributions 
from each subsample are also shown in Fig. El 

We use our simulation to estimate the backgrounds from NC and wrong-sign events. As 
shown in Fig. HJ the NC background is less than 2% for both neutrinos and antineutrinos. 
It increases with energy due to the contribution from high-inelasticity events in which the 
primary track is misidentified or perturbed by hits from the hadronic shower particles. The 
wrong-sign contamination is negligible in the neutrino sample but sizable in the antineutrino 
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sample, up to 5% at high energy. Wrong-sign background events at low energy come from 
U/j, CC events in the peak of the neutrino beam, while at higher energies, the through-going 
muons have increasingly larger bend radii, making their charge determination less certain. 



Flux Extraction 



The "low-i/" method [171 . Il8l | relies on the independence of the differential cross section, 



da 



; , with energy in the limit v — > 0. The differential dependence of the neutrino (antineu- 

prncc copfinn - 

dxdy 



trino) cross section, ^rjr-, on inelasticity, y = u/E, and the Bjorken scaling variable, x, can 



be written as 



<P<r"M G\ME ( 
dxdy 7i 



fx Mx v 2 A + (^)\i r v 2 i \ 



where Gf is the Fermi weak coupling constant, M is the proton mass and E is the incident 
neutrino energy. The plus sign in front of the xF 3 term is for neutrinos and the minus is 
for antineutrinos. The structure functions F 2 (x, Q 2 ), xF 3 (x, Q 2 ) and Rl(x, Q 2 ) depend on x 
and Q 2 , Rl is the ratio of the cross section for scattering from longitudinally to transversely 
polarized W-bosons. For quasi-elastic interactions the cross section can be written in this 
form with combinations of form factors replacing the structure functions. 

Integrating over x, the differential dependence on v can be written in the simplified form 

da v ' v / B v C v 2 \ 

A l + T-iF-TTHFJ • ( 2 ) 



dv \ AE A2E 2 / 

The coefficients A, B, and C depend on integrals over structure functions 



G 2 F M 



A = — £ — / F 2 (x)dx, (3) 

7T 
i2 



G 2 MP/ \ 
B = / [F 2 (x)TxFz(x)yx, 

g 2 mt 

C = B- / F 2 (x)Rdx, 

7T ./ 



where 
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The factor A is nearly the same for neutrino and antineutrino probes 1 , however, the mag- 
nitude of the coefficient B is larger for antineutrinos, where the xF% contribution is added, 
compared with the neutrino case where the term is subtracted. As discussed later, this 
makes the energy dependence correction needed in this method larger for the antineutrino 
flux shape. The C term, which depends on R^, is small. 

For small u/E, Eq. [2] shows that the differential cross section becomes independent of 
energy and is equal to the same constant, A, for neutrinos and antineutrinos. Multiplying 

both sides by the flux, (f>(E), and taking the limit u — > gives 

dN 

— = A*(E). (4) 

dv i/^o 

Therefore, the flux in a given energy bin can be approximated using the number of events 
at low v. 

We account for the small u/E and (u/E) 2 dependence resulting from a finite z/ in Eq. [2] 
using a "low-z/" correction 

a(u < u ,E — >• oo) 

that is calculated from our cross section model. The term a[y < uq, E) is the value of the 
integrated cross section below our chosen u cut at energy E, and o[y < u , E — >■ oo) is its 
value in the high energy limit. 

This correction is applied to our selected flux sample, F U ( U \E), consisting of the sub- 
set of Nqq\e) with v < vq that is subsequently corrected for acceptance, A^ U \E), and 
backgrounds, B^ u \e), 

S"<">K E) x Af \E) 

This yields the shape of the flux with energy. A normalization factor H u , determined using 
external data, must be applied to give the absolute flux, $> V ( U \E) = H v <\> v{jJ \E), as described 
in the next section. 

Our choice of uq trades statistical precision for modeling uncertainty in determining 
S v ^ u \uo, E). To improve the statistical precision, we increase uq with energy while keep- 
ing the ratio u/E and the resulting model dependence small. We set u =1 GeV for events 



1 For an isoscalar target with only u and d quarks F% = F%, assuming isospin symmetry. Including s 
quarks and CKM mixing gives a small difference term F% — F% — —\V% s (u v + d v ), where u v and d v 
are the valence quark distributions. We apply a correction to account for this term to the antineutrino 
normalization. 
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FIG. 5: The low-v correction, S u ^\uq, E), applied to the flux sample for neutrinos (left) and 
antineutrinos (right). The solid line shows S v ^(vq,E) for v <lGeV applied from 3-9 GeV (5- 
9 GeV for 17), dashed line for v <2GeV applied from 9-18 GeV, and the dotted line for v <5GeV 
applied above 18 GeV. 

with E v <9GeV, v Q =2 GeV for 9< E u <18GeV and u =5 GeV for E v >18GeV. Fig. 
shows the size of the \ow-u correction for neutrino and antineutrino samples. The correction 
for neutrinos is about 3% at 3 GeV and for antineutrino is about 20% at 5 GeV. 

The stronger inelasticity dependence of the antineutrino cross section results in the much 
larger correction for antineutrinos. In addition, antineutrino CC interactions have lower 
inelasticity on average, which causes a large overlap between the cross section and the flux 
samples. The overlap decreases with energy from 90% at 3 GeV to about 60% at 6 GeV for 
antineutrinos, whereas for neutrinos it is 60% at 3 GeV and below 30% above 6 GeV 44]. 
We therefore restrict our analysis to the region above 5 GeV for the antineutrino sample. 

The low-v correction introduces a model dependence and model uncertainty to the flux 
determination. We account for this uncertainty in the flux by varying the model parameters 
described in Sec. Ill CI and re-calculating the flux. The change in the correction when the 
model is varied is 1% or less because it is a fractional term in which the numerator and 
denominator are similarly affected. 

Fig. [6] and Table [Til show the extracted neutrino and antineutrino fluxes in the selected 
fiducial volume after normalization. The systematic uncertainties on the extracted flux are 
discussed in Sec. IIVI We correct the input flux model shown in Fig. [1] by reweighting the 
simulation with the ratio of the extracted flux to the original simulated flux. The resulting 
corrections to the initial simulated flux are consistent with those obtained by a different 
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E bin 


v Flux 


Error 


v Flux 


Error 


(GeV) 


Particles/GeV/m 2 /10 9 PoT 


3-4 


8 05xl0 4 


5.2xl0 3 






4-5 


3 06xl0 4 


2.4xl0 3 






5-7 


9 07xl0 3 


5.3 xlO 2 


2.80xl0 3 


330 


7-9 


5.18xl0 3 


3.5xl0 2 


2.32xl0 3 


170 


9-12 


3.21xl0 3 


2.2xl0 2 


1.32xl0 3 


85 


12-15 


1.94xl0 3 


l.OxlO 2 


6 89xl0 2 


42 


15-18 


1.09xl0 3 


65 


3 79xl0 2 


24 


18-22 


629 


37 


190 


14 


22-26 


348 


20 


86.3 


7.8 


26-30 


200 


13 


40.1 


3.9 


30-36 


119 


6.8 


19.3 


1.9 


36-42 


72.2 


3.9 


9.6 


0.9 


42-50 


51.6 


2.8 


4.9 


0.5 



TABLE II: Measured flux as a function of neutrino energy. Statistical, systematic and normaliza- 
tion uncertainties are included in the error estimate. 



technique |27j used for the MINOS oscillation analyses. 

Fig. [7] shows a comparison of the CC data sample and Monte Carlo simulation before 
and after flux reweighting is applied for the measured kinematic variables; the muon energy, 
hadronic shower energy, and muon track angle with respect to the beam direction. The 
agreement between Monte Carlo and data in all three distributions is significantly improved 
after the flux reweighting has been applied. 



D. Cross Section Extraction 



The energy dependence of the cross section is extracted by dividing the selected CC 
rate Tqq(E) by the measured flux & U ( U \E) in each energy bin. Explicitly including the 
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FIG. 6: Extracted flux as a function of neutrino energy. The error bars are too small to see on this 
log-scale plot, but are given in Table ITU 



corrections described above, the cross section is 
a$g(E) _ 1 r"£\E) _ 1 / (N" C P(E) -B&\E))/A$g(E) 



(7) 



E E&W{E) E \s»W(v ,E) xH-x (F"<P\E) - b£>\e)) / A£ p) {E) 
Because our data extends above 30 GeV we can take advantage of existing high en- 
ergy precision measurements [lisl, Q] to determine the flux normalization constant 
H u . H u is chosen so that our measured flux results in an average CC cross section 
value from 30-50 GeV that agrees with the world average on an isoscalar target (0.675 ± 
0.009) xlO~ 38 cm 2 /GeV in the same energy range. The uncertainty on the normalization 
constant from the statistical precision of our data and that from the world average cross sec- 
tion are added in quadrature and applied as an uncertainty on our measured cross section at 
all energies. We apply a normalization correction to the antineutrino sample to account for 
the small difference in i 7 ^ and F%. The correction, which is computed from our cross section 



model, is 1% for u <1 GeV, 2.6% for v <2 GeV and 3.8% for the v <5 GeV sample (44|. 

After extracting the flux with the low-i/ method, the cross section analysis is then repeated 
with the measured flux as input. This removes the effect of inaccuracies in the initial 
simulated flux on the acceptance corrections that are applied to both the flux and the cross 
section samples. The change in cross section between the final value and that extracted with 
the default simulated flux is less than 0.5% averaged over all data points. Since this effect 
is small we do not iterate the procedure further. 
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FIG. 7: Comparison of data and Monte Carlo simulation distributions for the kinematic variables 
Eft, Eh a d, and 9^ for neutrinos (left) and antineutrinos (right). The points show the data, the 
dashed line shows the nominal Monte Carlo model and the solid line shows the Monte Carlo 
model after applying flux reweighting. The ratio of the data to Monte Carlo is shown below each 
distribution. The error bars show the statistical and systematic uncertainties added in quadrature. 
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The cross section we report is that expected for an isoscalar target. The MINOS iron- 
scintillator detector has a 6.1% excess of neutrons over protons for which we correct using 
the NEUGEN3 cross section model 23]. The energy dependent corrections are about -2% for 
neutrinos and +2% for antineutrinos. We also apply corrections for radiative effects 45], 
which have an effect on the result of less than 1% at all energies. 

Appendix I VII provides the measured raw ratio of cross section to flux data samples where 
both numerator and denominator are corrected only for detector effects and backgrounds. 
As described further therein, this allows one to remove the model dependence in the cross 
section extraction and to use an alternative cross section model to correct the ratio and 
compute cross sections. 



IV. SYSTEMATIC UNCERTAINTIES 



Systematic uncertainties in the measurement from the following sources were considered: 
reconstructed muon and hadron energy scales, the effect of final state interactions on the 
measured energy, NC contamination, wrong-sign contamination (antineutrino sample only), 
our lack of perfect understanding of the detector and event reconstruction, and cross section 
modeling. Each systematic uncertainty is evaluated independently and propagated through 
the analysis, including a recalculation of the absolute normalization of the result. Many 
systematic effects cause changes that are similar in the cross section and flux samples and 
therefore partially cancel in the measured cross section. 

The largest uncertainty comes from knowledge of the absolute muon and hadronic energy 
scales discussed in Sec. Ill Bt The muon energy scale is more important for the flux sample 
than for the cross section sample because a larger fraction of the neutrino energy per event 
is carried by the muon in the former. Conversely, the hadronic energy scale is less important 
in the flux sample. 

Fig. [8] shows the effects of muon and hadronic energy scale uncertainties on the extracted 
cross section. These are evaluated by applying the one-sigma shift in each scale factor 
to the data, extracting a new flux, and determining a new cross section, including a new 
normalization to the world data between 30 to 50 GeV. The resulting curves in Fig. [8] 
represent the change from the baseline cross section measurement. These uncertainties peak 
in the low energy region and become small at high energies where the normalization is 
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pinned. 

The energy dependence of these uncertainties has a non-trivial shape because of the 
interplay of the shape of the flux spectrum and the method of normalizing to the world 
high energy measurements. In particular, applying a muon energy scale factor shifts the 
observed flux in one direction and causes inflection points in the shape of the extracted 
flux near 6 GeV and 14 GeV for neutrinos, which propagate (with some cancellation) to the 
cross section analysis as shown in Fig. [HJ A similar effect arises more directly in the cross 
section sample for the hadron energy scale. We show the +lcr systematic alone so that the 
shape distortion is clear. The -la distortion is the same shape but inverted in sign. The 
antineutrino analysis has fewer inflections because the spectrum is not peaked. 

As described in Sec. Ill C[ final state interactions affect the measured hadronic energy and 
are included in our Monte Carlo simulation. The uncertainty in their modeling contributes an 
effective hadronic energy scale uncertainty of 8% below hadronic energy of 1 GeV, decreasing 
to 4% above 5 GeV. Their effect on the cross section, shown in Fig. [BJ peaks at low energy 
and is fractionally larger for antineutrinos, which have a larger fraction of low hadron energy 
events due to their inelasticty distribution. 

The uncertainty from our knowledge of the NC contamination is obtained by varying 
the value of the minimum requirement, which selects the CC sample, from its nominal 
value of 1.5 GeV up to 2.0 GeV and down to 1.0 GeV. The resulting change in a/E is small 
and corresponds to a change of less than 1%, which we take to be the NC contamination 
uncertainty. 

To account for uncertainties in the acceptance correction, which arise from modeling of 
detector geometry, alignment, and magnetic field, we collected a dedicated data set with 
the detector magnetic field polarity reversed (set to focus positive charges) from its nominal 
running mode which focuses negative charges. In this data set, muon tracks pass through a 
different region of the detector and previously focused tracks bend away from the coil region. 
The flux and cross section are extracted and are compared with their values measured in 
the nominal mode. The flux sample relies more heavily on the muon track measurement 
than does the cross section sample and consequently larger effects are seen in the measured 
flux (on average 5% for neutrino and 2.5% for antineutrino flux). Smaller differences are 
expected for the antineutrino sample because more of the tracks exit the detector in the 
downstream region and do not pass near the difficult-to-model coil region in either mode. A 
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systematic uncertainty on the nominal polarity flux is taken to be half the difference between 
it and the flux extracted in the reversed polarity sample. This is added in quadrature with 
the other uncertainties. The cross section sample relies less heavily on the track momentum. 
The differences in measured cross section are at the level of 1%, which are neglected. 

Since the cross section model is used to apply a small energy dependent correction to 
the flux sample (see Eq. [5]), we take into account uncertainties in the model parameters 
described in Sec. Ill CI We account for uncertainties in the quasi-elastic and resonance 
contributions to the cross section by varying the axial mass parameters, M® E and M^ ES 
in our model by ±15%. The resulting effect on the cross section is less than 2%. In this 
measurement, we include an additional uncertainty in the DIS component of the model to 
account for contributions to the v dependence of the cross section that could affect the flux 
extraction. To quantify the resulting uncertainty, we vary each parameter in the model 34] 



and study the change of the reduced x 2 °f the fit to the charged lepton data from which 
they were originally determined. We take the shift that corresponds to a one unit shift in 
fit x 2 as the uncertainty for each parameter. The values and the associated uncertainties of 
A ht , B ht , C vlu and C v2u are determined to be 0.538±0.134, 0.305±0.076, 0.291±0.087 and 
0.189±0.076, respectively. The other parameters (C v id, C v2 d, C su and C s d) have a negligible 
effect on the analysis. The effects of these uncertainties on the flux measurement and the 
acceptance correction are propagated to an uncertainty on the extracted cross section. The 
contributions from each parameter shift are added linearly to form the total DIS model 
uncertainty which is 2% below 8 GeV for both neutrino and antineutrino cross sections, and 
is negligible above this energy. 

The contamination from wrong-sign events is significant only for the antineutrino sample. 
To evaluate the uncertainty from this source, we recompute the cross section assuming no 
wrong-sign contamination and two times as much wrong-sign contamination. The resulting 
uncertainty is negligible below 15 GeV but is about 4% at the highest energies. 

The systematic, statistical, and total uncertainties, for the neutrino and antineutrino 
cross sections are summarized in Fig. [9j 

In the cross section ratio significant additional cancellation of uncertainties occurs. At the 
lower energies the uncertainties are about half those shown in Fig. [HJ At higher energies they 
are one to two percent each, except for the wrong-sign contamination which is significant 
only for the antineutrino sample. 
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FIG. 8: Effect of the energy scale uncertainty on the neutrino (left) and antineutrino (right) 
extracted cross section. The curves give the shape distortion due to a one-sided one-sigma error. 
The solid line shows the effect of increasing the muon energy scale by 2% for stopping muons and 
4% for exiting muons. The dotted line shows the effect of increasing the hadronic energy scale by 
5.6%, and the dashed line shows the effect of shifting the final state interaction model. The plot 
on the right also shows the effect of wrong-sign contamination uncertainty (large dashes) on the 
antineutrino cross section. 




FIG. 9: Summary of the statistical, total systematic, and total uncertainty for the neutrino (left) 
and antineutrino (right) extracted cross section. 
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V. RESULTS 



Fig. [10] shows the extracted energy dependence of the total cross section divided by 
energy (a/E) for v^N CC and for v^N CC interactions on an isoscalar target. The cross 
section values are assigned to the average energy in the bin. Both cross sections approach a 
linear energy dependence for energies above 20 GeV. For neutrinos a/E drops with increasing 
energy in the lower energy region. At 3 GeV the quasi-elastic cross section is still expected 
to be appreciable (~15%). Its contribution to a/E falls rapidly with increasing energy as 
inelastic processes (resonance production and DIS) turn on. For antineutrinos the measured 
a/E rises gradually in the region 5-20 GeV to its asymptotic high energy value. In this 
case the falling fractional contribution of the quasi-elastic cross section is offset by the more 
gradual turn on of the DIS process, which is expected due to its strong dependence on 
the antiquark component which rises slowly with increasing Q 2 . Table [TTTJ summarizes the 
neutrino and antineutrino cross section results. 

Fig. [UJ shows MINOS neutrino and antineutrino results compared to the results from 
other experiments. The MINOS neutrino cross section agrees with previous measurements 
from CRS |g|, SKAT Q, IHEP-JINR 12], and GGM-PS 10 1 , but these experiments have 



significantly larger uncertainties. Our neutrino cross section is in good agreement and has 
comparable precision with the recent NOMAD measurement [15] . Our result is systematics- 
limited in the region below 15 GeV where the largest uncertainties comes from knowledge of 
the absolute muon and hadronic energy scales, whereas in NOMAD the flux determination 
dominates the uncertainty. The MINOS antineutrino cross section result is in good agree- 
ment with the sparse data available at lower energies and has much smaller uncertainty in 
the 10 to 30 GeV region. 

Fig. [12] and Table IIVI show the ratio of the P^N CC to i^N CC inclusive cross section 
as a function of energy. Because of cancellation of many of the systematic uncertainties the 
MINOS result is statistics-limited above 10 GeV. The cross section ratio appears to gradually 
approach its asymptotic scaling value of 0.504±0.003, defined by the world average from 30- 
200 GeV calculated from previous experiments [4]- The MINOS average ratio measured from 
30-50 GeV of 0.489±0.012 is in good agreement with the asymptotic value. At 10 GeV the 
measured ratio is 14% below the asymptotic value with 6.6a significance and at 24 GeV the 
measurement lies 7% below with 2.4<r significance. Our precise data show a slower approach 
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FIG. 10: Neutrino (left) and antineutrino (right) charged-current inclusive cross section per nucleon 
divided by energy for an isoscalar iron target. The black error bars show the statistical uncertainty 
and the shaded boxes show the statistical and systematic uncertainties added in quadrature. The 
dotted band shows the uncertainty on the normalization (about 1.5%). The solid black line shows 



the work 



re world average neutrino cross section value of 0.675 x 10 38 cm 2 /GeV from 30 to 50 GeV 



D- 



and the dashed black line shows this value extrapolated to lower energies. The neutrino 
cross section above 30 GeV is normalized using this world average value and the same normalization 
constant is then applied to the antineutrinos. The solid black line on the antineutrino cross section 
plot shows a world average antineutrino cross section value of 0.329 x 10 -38 cm 2 /GeV from 30 to 



50 GeV 
normalization. 



Ill ] . This value is shown for comparison and is not used for antineutrino sample 



to scaling behavior than has been previously claimed by low energy measurements 2 , which 



found their data to be consistent with scaling in the few GeV range 

Fig. [12] also shows the cross section ratio compared with the few other existing measure- 
ments. The MINOS data uniquely spans the 10-30 GeV region. It overlaps the precise high 



energy measurements 



13, 



2l| as well as the Gargamelle low energy measurement 22| which 



has precision of only about 20%. 

The total neutrino and antineutrino cross section in the Quark Parton Model, which 
describes neutrino scattering at high energy, can be written as 



a(uN) = ^^(Q+ 1 -Q), 

7T O 



(8) 



2 These are based on fits to the measured neutrino cross section with energy. 
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16.41 


0.304 


0.007 


0.007 


0.006 


0.012 


18-22 


19.87 


0.689 


0.006 


0.009 


0.012 


0.016 


19.82 


0.316 


0.006 


0.011 


0.005 


0.013 


22-26 


23.88 


0.683 


0.008 


0.005 


0.012 


0.015 


23.82 


0.320 


0.009 


0.004 


0.005 


0.011 


26-30 


27.89 


0.686 


0.010 


0.004 


0.012 


0.016 


27.84 


0.332 


0.012 


0.005 


0.006 


0.015 


30-36 


32.81 


0.675 


0.010 


0.002 


0.011 


0.016 


32.72 


0.325 


0.014 


0.006 


0.005 


0.016 


36-42 


38.87 


0.675 


0.013 


0.005 


0.011 


0.018 


38.74 


0.352 


0.021 


0.011 


0.006 


0.024 


42-50 


45.77 


0.676 


0.014 


0.004 


0.011 


0.019 


45.61 


0.324 


0.023 


0.013 


0.005 


0.027 



TABLE III: Summary of neutrino (left) and antineutrino (right) cross section results. The second 
column for each species shows the average energy in each bin. The uncertainties shown in columns 
4-7 for each species are the statistical, systematic, normalization and total contributions, respec- 
tively. The total uncertainty is obtained by summing the statistical, systematic, and normalization 
uncertainties in quadrature. 

and 

a(VN) = ^^(Q + 1 -Q) (9) 

71 6 

with the ratio of the two given by 

r _ g _ 1 + nm 
a{vN) 3 + Q/Q ' 1 ' 

where Q = J x[u(x) + d(x)]dx and Q = J x\u(x) +d(x)]dx. Here u(x) (u{x)) and d(x) (d(x)) 
are the parton distribution functions (PDFs) for the up and down quarks (antiquarks) in 
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FIG. 11: MINOS neutrino and antineutrino charged-current inclusive cross section compared with 



other experimental results 



8-13 



15, 



17 



19 21 



46 



47J. The error bars show the statistical, 



systematic, and normalization uncertainties added in quadrature. The solid black line shows the 
average world cross section in the 30 to 50 GeV region for the neutrino (0.675 x 10~ 38 cm 2 /GeV) 
and the antineutrino (0.329 x 10~ 38 cm 2 /GeV). The dashed line shows these high energy values 
extrapolated to lower energies. 

the nucleoli, respectively. In the limit of large Q 2 , Q 2 » M 2 , the PDFs depend only 
on x and are independent of Q 2 . In this limit the QPM predicts scaling behavior, i.e., 
a linear dependence of the cross sections with energy. In the low energy (low Q 2 ) limit, 
scaling violations occur and the QPM breaks down. Scattering off the entire nucleon (quasi- 
elastic scattering) and resonance production, where the nucleon is excited and decays to low 
multiplicity final states, must also be considered to account for the energy dependence of 
the cross section. 

The ratio r is constant with energy in the QPM and depends only on the integrated quark 
and antiquark distributions in the high-Q 2 limit. Eq. |9] indicates antiquarks are relatively 
more important in the antineutrino scattering case, r approaches the limiting value of 1/3 
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FIG. 12: (left) Ratio of antineutrino-nucleon to neutrino-nucleon cross section as a function of 
energy. Black error bars show the statistical uncertainty and shaded boxes show the total uncer- 
tainty with statistical and systematic uncertainties added in quadrature. The solid black line at 
0.504±0.003 is drawn at the average value obtained from previous measurements over the energy 



range 30-200 GeV 



J], (right) Comparison of measured r with other measurements for E <100 GeV. 



The MINOS result spans the intermediate energy r ang e and overlaps with the low energy data 



as well as with precise high energy measurements 



13, 



2l( 



22] 



if antiquarks are not present in the nucleon. High energy measurements of r can be used to 
measure the fraction of momentum carried by antiquarks in the nucleon. 

In order to interpret our measurement of r in the context of the QPM the quasi-elastic 
contribution is removed from the measured value by denning r ine i = o v inel j a v ind to be the 
cross section ratio for the purely inelastic contribution to the cross section. To compute 



r ine i, the NEUGEN3 cross section model 



23] is used to remove the fractional quasi-elastic 



contribution 3 . Table HVl gives the measured ratio r and the inelastic fraction r ine i along with 
their experimental uncertainties. The similarly slow increase of ri ne i with energy shows that 
the decrease in the quasi-elastic contributions alone has only a small effect on the observed 
shape. 

Eq. CEUlcan be rearranged to give the fraction of total quark momentum in the nucleon that 
is carried by antiquarks, = ~ ?^~n ■ This fraction as a function of energy is also given 



3 The neugen3 model, as described earlier, uses a value of M® E = 0.99. We provide the raw measured 
r so that one can use other models to compute the inelastic fraction. This will be especially useful as 
knowledge of M® E improves. For reference, increasing M® E by 0.15 decreases the inelastic fraction by 
less than 1% at 5.9 GeV, which is small compared with the experimental uncertainty. 
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Energy (GeV) 


r 


Stat. err. 


Total err. 






Q 
Q+Q 


A( $ ) 

^q+q' 


5.9 


0.444 


0.007 


0.032 


0.407 


0.029 


0.079 


0.030 


8.0 


0.417 


0.008 


0.016 


0.389 


0.016 


0.060 


0.016 


10.5 


0.433 


0.006 


0.010 


0.410 


0.010 


0.081 


0.010 


13.4 


0.454 


0.009 


0.010 


0.435 


0.010 


0.106 


0.010 


16.4 


0.430 


0.011 


0.012 


0.415 


0.012 


0.086 


0.012 


19.9 


0.457 


0.010 


0.015 


0.444 


0.014 


0.115 


0.014 


23.9 


0.467 


0.014 


0.015 


0.455 


0.015 


0.126 


0.014 


27.9 


0.482 


0.019 


0.022 


0.472 


0.021 


0.142 


0.019 


32.8 


0.480 


0.021 


0.023 


0.472 


0.023 


0.141 


0.021 


38.9 


0.520 


0.032 


0.037 


0.512 


0.037 


0.177 


0.032 


45.8 


0.477 


0.036 


0.041 


0.471 


0.040 


0.140 


0.037 



TABLE IV: The measured cross section ratio r at the bin average energy along with statistical and 



total uncertainties are given. To compute ri ne i, the neugen3 



231 ] cross section model is used to 



remove the fractional quasi-elastic contribution. The quark parton model is used to estimate the 
fraction of the total quark momentum that is carried by antiquarks, ■ Uncertainties computed 
for Ari ne i and the antiquark fraction do not include any model uncertainty contributions. 



in Table IIVI As neutrino energy decreases, one moves increasingly away from the domain 
of validity of this expression, which is derived in the DIS region. Target mass corrections 
as well as higher-twist terms become more important, especially at high-x. However, the 
high-x region contributes little to the Q and Q integrals. In the approximation that the 
contributions from these effects are small, our results are consistent with a non-zero antiquark 
content in the nucleon at our lowest energy, 5.9 GeV {(Q 2 ) = 1.4 GeV 2 ) and a gradual 
increase of the antiquark fraction with energy. In order to accurately extract the antiquark 
fraction from our data, a full higher order QCD model that incorporates these effects is 
required. 
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VI. CONCLUSION 



We have measured the charged- current neutrino-nucleus inclusive cross section in the 
energy range 3-50 GeV with a precision of 2-8% and the antineutrino-nucleus cross section 
from 5-30 GeV with a precision in the range 3-9%. The flux was determined by using a 
subsample of low-hadronic-energy events to measure the flux shape and the world average 
cross section above 30 GeV for normalization. This method was previously used at higher 
energies 17|, |21[ and here we have extended it down to 3 GeV. While the measurements are 



systematics-dominated, the overall systematic uncertainty benefits from partial cancellation 
in detector related systematic uncertainties that arise from measuring the flux and the CC 
event rate in the same detector. Both measurements impact the precision of total cross 
section measurements in the less than 30 GeV range. 

Our measurement of the antineutrino to neutrino cross section ratio is the most precise in 



the less than 30 GeV range, where only one previous measurement has been performed [22]. 
The measured rise of the cross section ratio with energy is consistent with an expected 
slow rise in the antineutrino inelastic cross section with the increase in number of sea-quark 
degrees of freedom for increasing Q 2 . 

The measurement presented here can be used to tune neutrino and antineutrino cross 
section models which benefit ongoing and future neutrino oscillation measurements. 
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Appendix: Uncorrected Data Sample 

Our measurement has cross section model dependence which arises from the correction for 
the minimum muon energy requirement >1.5 GeV in the cross section sample (Nqq\e)) 
and for the small energy dependence in the low-v flux sample (see Eq. [2]). Here we provide 
the raw data cross section to flux sample ratio R u ^ corrected only for detector effects and 
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backgrounds. This will allow the reader to use an alternative cross section model and our 
data to compute neutrino and antineutrino cross sections. 

Table |V] gives the values of R v (R u for antineutrinos) where both numerator and de- 
nominator have been corrected for detector effects and backgrounds using our geant3- 
based detector simulation. The ratio has not been corrected for the effect of the kinematic 
E a >1.5GeV cut, which affects only the numerator. The unnormalized neutrino cross sec- 
tion o v unnorm {E) can be computed from R U (E) by applying two corrections, 

<nnorm(E) = K\E) X S> , E) X R»(E), (11) 

where K V [E) = N(E)/N(E,E fl >1.5GeV) is the ratio of total cross section events for all 
muon energies to that with muon energies larger than 1.5 GeV in each energy bin, and 
S u (u , E) is defined in Eq. [5j (where vq=1,2 or 5 GeV). 

The cross section is normalized using the values of R v for different v§ cut samples provided 
in Table IVII To improve statistical precision of the flux sample three different values of the 
uq cut were used in the analysis; uq <lGeV applies for the neutrino energies E <9GeV, 
vq <2 GeV for 9< E <18 GeV, and u <5 GeV for E >18 GeV. We define the normalization 
constant Norm(u ) for each v sample as 

Korm( E ) = Norm(u ) x a v unnorm (E). (12) 

Norm{vo) is obtained by first computing unnormalized cross sections in the range 30< 
E v <50GeV using R v from Table [VTl and Eq. [TTJ The weighted average of a v unnorm {Ei) / (E v ) 
in this energy range is computed for each Uq sample using the statistical errors on R v also 
given in Table IVII Normiv^) is then obtained by scaling this to the world average value 
0.675xl(T 38 cm 2 /GeV. 

The same overall normalization constants are used to obtain the antineutrino cross section 

al^jE) = [Norm(u ) x G corr (u )} x K V (E) x S> , E) x R V {E). (13) 

The additional normalization factor G corr {vo), is used to account for a small difference in 
neutrino and antineutrino F 2 structure functions (see Sec. IIII CJT) . G corr (z/ ) can be computed 
from the ratio of asymptotic values of antineutrino to neutrino low-v cross sections, 

v u; o v (v < u ,E oo) K J 
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Neutrino 


Antineutrino 


E bin (GeV) 


(E u ) (GeV) 


R" 


{Eh) (GeV) 


R" 


3-4 


3.48 


1.72 






4-5 


4.45 


2.35 






5-7 


5.89 


3.31 


6.07 


2.22 


7-9 


7.97 


4.93 


7.99 


2.78 


9-12 


10.5 


3.71 


10.4 


2.10 


12-15 


13.4 


4.87 


13.4 


2.69 


15-18 


16.4 


6.22 


16.4 


3.11 


18-22 


19.9 


3.32 


19.8 


1.90 


22-26 


23.9 


3.98 


23.8 


2.22 


26-30 


27.9 


4.68 


27.8 


2.60 


30-36 


32.8 


5.45 


32.7 


2.94 


36-42 


38.9 


6.49 


38.7 


3.75 


42-50 


45.8 


7.70 


45.6 


4.04 



TABLE V: Ratio of cross section to flux sample where both numerator and denominator have been 
corrected for detector effects and backgrounds. The text describes how to use this data and a cross 
section model to compute neutrino and antineutrino cross sections. 

for the three different z/ values. 

The uncertainty in the new measured cross section should be estimated using the frac- 
tional cross section uncertainty (syst©stat) given in Table IHIl which properly takes into 
account cancellations in several systematic uncertainties in the cross section and flux sam- 
ples. 
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